TUT4/7-mediated uridylation of a coronavirus subgenomic RNAs delays viral replication.

Coronaviruses are positive-strand RNA viruses with 3' polyadenylated genomes and subgenomic transcripts. The lengths of the viral poly(A) tails change during infection by mechanisms that remain poorly understood. Here, we use a splint-ligation method to measure the poly(A) tail length and poly(A) terminal uridylation and guanylation of the mouse hepatitis virus (MHV) RNAs. Upon infection of 17-CL1 cells with MHV, a member of the Betacoronavirus genus, we observe two populations of terminally uridylated viral transcripts, one with poly(A) tails ~44 nucleotides long and the other with poly(A) tails shorter than ~22 nucleotides. The mammalian terminal uridylyl-transferase 4 (TUT4) and terminal uridylyl-transferase 7 (TUT7), referred to as TUT4/7, add non-templated uracils to the 3'-end of endogenous transcripts with poly(A) tails shorter than ~30 nucleotides to trigger transcript decay. Here we find that depletion of the host TUT4/7 results in an increased replication capacity of the MHV virus. At late stages of infection, the population of uridylated subgenomic RNAs with tails shorter than ~22 nucleotides is reduced in the absence of TUT4/7 while the viral RNA load increases. Our findings indicate that TUT4/7 uridylation marks the MHV subgenomic RNAs for decay and delays viral replication.

Prophylactic intranasal administration of lipid nanoparticle formulated siRNAs reduce SARS-CoV-2 and RSV lung infection.

RNA interference (RNAi) is an emerging and promising therapy for a wide range of respiratory viral infections. This highly specific suppression can be achieved by the introduction of short-interfering RNA (siRNA) into mammalian systems, resulting in the effective reduction of viral load. Unfortunately, this has been hindered by the lack of a good delivery system, especially via the intranasal (IN) route. Here, we have developed an IN siRNA encapsulated lipid nanoparticle (LNP) in vivo delivery system that is highly efficient at targeting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and respiratory syncytial virus (RSV) lung infection in vivo. Importantly, IN siRNA delivery without the aid of LNPs abolishes anti-SARS-CoV-2 activity in vivo. Our approach using LNPs as the delivery vehicle overcomes the significant barriers seen with IN delivery of siRNA therapeutics and is a significant advancement in our ability to delivery siRNAs. The study presented here demonstrates an attractive alternate delivery strategy for the prophylactic treatment of both future and emerging respiratory viral diseases.

Regulation of ribonucleoprotein condensates by RNase L during viral infection.

In response to viral infection, mammalian cells activate several innate immune pathways to antagonize viral gene expression. Upon recognition of viral double-stranded RNA, protein kinase R (PKR) phosphorylates the alpha subunit of eukaryotic initiation factor 2 (eIF2α) on serine 51. This inhibits canonical translation initiation, which broadly antagonizes viral protein synthesis. It also promotes the assembly of cytoplasmic ribonucleoprotein complexes termed stress granules (SGs). SGs are widely thought to promote cell survival and antiviral signaling. However, co-activation of the OAS/RNase L antiviral pathway inhibits the assembly of SGs and promotes the assembly of an alternative ribonucleoprotein complex termed an RNase L-dependent body (RLB). The formation of RLBs has been observed in response to double-stranded RNA, dengue virus infection, or SARS-CoV-2 infection. Herein, we review the distinct biogenesis pathways and properties of SGs and RLBs, and we provide perspective on their potential functions during the antiviral response. This article is categorized under: RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Turnover and Surveillance > Regulation of RNA Stability RNA Export and Localization > RNA Localization.

Remdesivir does not affect mitochondrial DNA copy number or deletion mutation frequency in aged male rats: A short report.

Remdesivir is a leading therapy in patients with moderate to severe coronavirus 2 (SARS-CoV-2) infection; the majority of whom are older individuals. Remdesivir is a nucleoside analog that incorporates into nascent viral RNA, inhibiting RNA-directed RNA polymerases, including that of SARS-CoV-2. Less is known about remdesivir's effects on mitochondria, particularly in older adults where mitochondria are known to be dysfunctional. Furthermore, its effect on age-induced mitochondrial mutations and copy number has not been previously studied. We hypothesized that remdesivir adversely affects mtDNA copy number and deletion mutation frequency in aged rodents. To test this hypothesis, 30-month-old male F333BNF1 rats were treated with remdesivir for three months. To determine if remdesivir adversely affects mtDNA, we measured copy number and mtDNA deletion frequency in rat hearts, kidneys, and skeletal muscles using digital PCR. We found no effects from three months of remdesivir treatment on mtDNA copy number or deletion mutation frequency in 33-month-old rats. These data support the notion that remdesivir does not compromise mtDNA quality or quantity at old age in mammals. Future work should focus on examining additional tissues such as brain and liver, and extend testing to human clinical samples.

Evaluation of the mutagenic effects of Molnupiravir and N4-hydroxycytidine in bacterial and mammalian cells by HiFi sequencing.

Molnupiravir (MOV) is used to treat COVID-19. In cells, MOV is converted to the ribonucleoside analog N4-hydroxycytidine (NHC) and incorporated into the SARS-CoV-2 RNA genome during its replication, resulting in RNA mutations. The widespread accumulation of such mutations inhibits SARS-CoV-2 propagation. Although safety assessments by many regulatory agencies across the world have concluded that the genotoxic risks associated with the clinical use of MOV are low, concerns remain that it could induce DNA mutations in patients, particularly because numerous in vitro studies have shown that NHC is a DNA mutagen. In this study, we used HiFi sequencing, a technique that can detect ultralow-frequency substitution mutations in whole genomes, to evaluate the mutagenic effects of MOV in E. coli and of MOV and NHC in mouse lymphoma L5178Y cells and human lymphoblastoid TK6 cells. In all models, exposure to these compounds increased genome-wide mutation frequencies in a dose-dependent manner, and these increases were mainly composed of A:T → G:C transitions. The NHC exposure concentrations used for mammalian cells were comparable to those observed in the plasma of humans who received clinical doses of MOV.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans.

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV's roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species' genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates' germlines after the segregation of Gibbons' species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold-fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells' populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein-protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer's, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells' differentiation programs.

Gene networks under circadian control exhibit diurnal organization in primate organs.

Mammalian organs are individually controlled by autonomous circadian clocks. At the molecular level, this process is defined by the cyclical co-expression of both core transcription factors and their downstream targets across time. While interactions between these molecular clocks are necessary for proper homeostasis, these features remain undefined. Here, we utilize integrative analysis of a baboon diurnal transcriptome atlas to characterize the properties of gene networks under circadian control. We found that 53.4% (8120) of baboon genes are oscillating body-wide. Additionally, two basic network modes were observed at the systems level: daytime and nighttime mode. Daytime networks were enriched for genes involved in metabolism, while nighttime networks were enriched for genes associated with growth and cellular signaling. A substantial number of diseases only form significant disease modules at either daytime or nighttime. In addition, a majority of SARS-CoV-2-related genes and modules are rhythmically expressed, which have significant network proximities with circadian regulators. Our data suggest that synchronization amongst circadian gene networks is necessary for proper homeostatic functions and circadian regulators have close interactions with SARS-CoV-2 infection.

When 3D genome technology meets viral infection, including SARS-CoV-2.

Mammalian chromosomes undergo varying degrees of compression to form three-dimensional genome structures. These three-dimensional structures undergo dynamic and precise chromatin interactions to achieve precise spatial and temporal regulation of gene expression. Most eukaryotic DNA viruses can invade their genomes into the nucleus. However, it is still poorly understood how the viral genome is precisely positioned after entering the host cell nucleus to find the most suitable location and whether it can specifically interact with the host genome to hijack the host transcriptional factories or even integrate into the host genome to complete its transcription and replication rapidly. Chromosome conformation capture technology can reveal long-range chromatin interactions between different chromosomal sites in the nucleus, potentially providing a reference for viral DNA-host chromatin interactions. This review summarized the research progress on the three-dimensional interaction between virus and host genome and the impact of virus integration into the host genome on gene transcription regulation, aiming to provide new insights into chromatin interaction and viral gene transcription regulation, laying the foundation for the treatment of infectious diseases.

Evolution of RNA sensing receptors in birds.

Birds are important hosts for many RNA viruses, including influenza A virus, Newcastle disease virus, West Nile virus and coronaviruses. Innate defense against RNA viruses in birds involves detection of viral RNA by pattern recognition receptors. Several receptors of different classes are involved, such as endosomal toll-like receptors and cytoplasmic retinoic acid-inducible gene I-like receptors, and their downstream adaptor proteins. The function of these receptors and their antagonism by viruses is well established in mammals; however, this has received less attention in birds. These receptors have been characterized in a few bird species, and the completion of avian genomes will permit study of their evolution. For each receptor, functional work has established ligand specificity and activation by viral infection. Engagement of adaptors, regulation by modulators and the supramolecular organization of proteins required for activation are incompletely understood in both mammals and birds. These receptors bind conserved nucleic acid agonists such as single- or double-stranded RNA and generally show purifying selection, particularly the ligand binding regions. However, in birds, these receptors and adaptors differ between species, and between individuals, suggesting that they are under selection for diversification over time. Avian receptors and signalling pathways, like their mammalian counterparts, are targets for antagonism by a variety of viruses, intent on escape from innate immune responses.

Poly(I:C) exposure during in vitro fertilization disrupts first cleavage of mouse embryos and subsequent blastocyst development.

The reproductive system can be infected by a variety of double-stranded RNA viruses, which disrupt ovary function and pregnancy. However, whether viral infection directly affects early embryonic development remains unknown. Here we show that Poly(I:C), which mimics a double-stranded RNA virus, significantly impaired mouse early embryonic development in vitro, and up-regulated TLR3 and IFNα at the two cells embryo stage. Further studies indicated that Poly(I:C)-treatment caused DNA damage and abnormal spindle morphology at the first cleavage. Moreover, CDX2 and SOX2 expression was decreased while blastocyst cell apoptosis was increased. Altogether, Poly(I:C) decreased the rate of successful in vitro fertilization via DNA damage and abnormal spindle morphology at the first cleavage and inhibited early embryonic development by inducing immune response and promoting blastocyst cell apoptosis. This study provides an implication for exploring the causes of reproductive disorders in mammals and humans caused by infection of double-stranded RNA virus.

Mammalian viral suppressors of RNA interference.

The antiviral defense directed by the RNAi pathway employs distinct specificity and effector mechanisms compared with other immune responses. The specificity of antiviral RNAi is programmed by siRNAs processed from virus-derived double-stranded RNA by Dicer endonuclease. Argonaute-containing RNA-induced silencing complex loaded with the viral siRNAs acts as the effector to mediate specific virus clearance by RNAi. Recent studies have provided evidence for the production and antiviral function of virus-derived siRNAs in both undifferentiated and differentiated mammalian cells infected with a range of RNA viruses when the cognate virus-encoded suppressor of RNAi (VSR) is rendered nonfunctional. In this review, we discuss the function, mechanism, and evolutionary origin of the validated mammalian VSRs and cell culture assays for their identification.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

Discovery and Use of Long dsRNA Mediated RNA Interference to Stimulate Antiviral Protection in Interferon Competent Mammalian Cells.

In invertebrate cells, RNA interference (RNAi) acts as a powerful immune defense that stimulates viral gene knockdown thereby preventing infection. With this pathway, virally produced long dsRNA (dsRNA) is cleaved into short interfering RNA (siRNA) by Dicer and loaded into the RNA-induced silencing complex (RISC) which can then destroy/disrupt complementary viral mRNA sequences. Comparatively, in mammalian cells it is believed that the type I interferon (IFN) pathway is the cornerstone of the innate antiviral response. In these cells, dsRNA acts as a potent inducer of the IFN system, which is dependent on dsRNA length, but not sequence, to stimulate an antiviral state. Although the cellular machinery for RNAi is intact and functioning in mammalian cells, its role to trigger an antiviral response using long dsRNA (dsRNAi) remains controversial. Here we show that dsRNAi is not only functional but has a significant antiviral effect in IFN competent mammalian cells. We found that pre-soaking mammalian cells with concentrations of sequence specific dsRNA too low to induce IFN production could significantly inhibit vesicular stomatitis virus expressing green fluorescent protein (VSV-GFP), and the human coronaviruses (CoV) HCoV-229E and SARS-CoV-2 replication. This phenomenon was shown to be dependent on dsRNA length, was comparable in effect to transfected siRNAs, and could knockdown multiple sequences at once. Additionally, knockout cell lines revealed that functional Dicer was required for viral inhibition, revealing that the RNAi pathway was indeed responsible. These results provide the first evidence that soaking with gene-specific long dsRNA can generate viral knockdown in mammalian cells. We believe that this novel discovery provides an explanation as to why the mammalian lineage retained its RNAi machinery and why vertebrate viruses have evolved methods to suppress RNAi. Furthermore, demonstrating RNAi below the threshold of IFN induction has uses as a novel therapeutic platform, both antiviral and gene targeting in nature.

First report and genome sequencing of SARS-CoV-2 in a cat (Felis catus) in Colombia.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a virus of zoonotic origin that can bind to ACE2 receptors on the cells of many wild and domestic mammals. Studies have shown that the virus can circulate among animals mutate, lead to animal-to-human zoonotic jump, and further onward spread between humans. Infection in pets is unusual, and there are few human-to-pet transmission reports worldwide. OBJECTIVE: To describe the SARS-CoV-2 infection in a domestic animal in Córdoba, Colombian Caribbean region. METHODS: A cross-sectional molecular surveillance study was carried out, oral and rectal swabs were taken from cats and dogs living with people diagnosed with coronavirus disease 2019 (COVID-19). RESULTS: SARS-CoV-2 was found in a cat living with a person with COVID-19. Genome sequencing showed that the B.1.111 lineage caused the infection in the cat. The owner's sample could not be sequenced. The lineage is predominant in Colombia, and this variant is characterised by the presence of the D614D and Q57H mutation. CONCLUSION: The present work is the first report of an infected cat with SARS-CoV-2 with whole-genome sequencing in Colombia. It highlights the importance of detecting SARS-CoV-2 mutations that could promote the transmissibility of this new coronavirus. There is still a significant information gap on human-to-cat-to-human infection; we encourage self-isolation measures between COVID-19 patients and companion animals. The findings of this study give a preliminary view of the current panorama of SARS-CoV-2 infection in animals in Colombia.

Adenovirus entry: Stability, uncoating, and nuclear import.

Adenoviruses (AdVs) are widespread in vertebrates. They infect the respiratory and gastrointestinal tracts, the eyes, heart, liver, and kidney, and are lethal to immunosuppressed people. Mastadenoviruses infecting mammals comprise several hundred different types, and many specifically infect humans. Human adenoviruses are the most widely used vectors in clinical applications, including cancer treatment and COVID-19 vaccination. AdV vectors are physically and genetically stable and generally safe in humans. The particles have an icosahedral coat and a nucleoprotein core with a DNA genome. We describe the concept of AdV cell entry and highlight recent advances in cytoplasmic transport, uncoating, and nuclear import of the viral DNA. We highlight a recently discovered "linchpin" function of the virion protein V ensuring cytoplasmic particle stability, which is relaxed at the nuclear pore complex by cues from the E3 ubiquitin ligase Mind bomb 1 (MIB1) and the proteasome triggering disruption. Capsid disruption by kinesin motor proteins and microtubules exposes the linchpin and renders protein V a target for MIB1 ubiquitination, which dissociates V from viral DNA and enhances DNA nuclear import. These advances uncover mechanisms controlling capsid stability and premature uncoating and provide insight into nuclear transport of nucleic acids.

Developing mRNA for Therapy.

Messenger RNA was discovered in 1961 and it took 60 years until the first mRNA became FDA-approved product in the form of COVID-19 mRNA vaccine. During those years a lot of progress has been made by hundreds of scientists. It was 1978 when the first-time isolated mRNA delivered into mammalian cells produced the encoded protein. In vitro transcription introduced in 1984 made it possible to generate any desired mRNA from the encoding plasmid using phage RNA polymerases. In the early 90s mRNA was used for therapy as well as vaccine against infectious diseases and cancer. Inflammatory nature of the mRNAs limited their in vivo use. Replacing uridine with pseudouridine made the mRNA non-immunogenic, more stable and highly translatable. Delivery of the lipid nanoparticle-formulated nucleoside-modified mRNA encoding viral antigens became a platform for effective vaccine. Labile nature of the mRNA is ideal for transient production of the viral antigen, to generate effective antibody and cellular immune response. The mRNA platform is revolutionizing the delivery of effective and safe vaccines, therapeutics and gene therapies.

Sindbis Macrodomain Poly-ADP-Ribose Hydrolase Activity Is Important for Viral RNA Synthesis.

ADP-ribosylation is a highly dynamic posttranslational modification frequently studied in stress response pathways with recent attention given to its role in response to viral infection. Notably, the alphaviruses encode catalytically active macrodomains capable of ADP-ribosylhydrolase (ARH) activities, implying a role in remodeling the cellular ADP-ribosylome. This report decouples mono- and poly-ARH contributions to macrodomain function using a newly engineered Sindbis virus (SINV) mutant with attenuated poly-ARH activity. Our findings indicate that viral poly-ARH activity is uniquely required for high titer replication in mammalian systems. Despite translating incoming genomic RNA as efficiently as WT virus, mutant viruses have a reduced capacity to establish productive infection, offering a more complete understanding of the kinetics and role of the alphavirus macrodomain with important implications for broader ADP-ribosyltransferase biology. IMPORTANCE Viral macrodomains have drawn attention in recent years due to their high degree of conservation in several virus families (e.g., coronaviruses and alphaviruses) and their potential druggability. These domains erase mono- or poly-ADP-ribose, posttranslational modifications written by host poly-ADP-ribose polymerase (PARP) proteins, from undetermined host or viral proteins to enhance replication. Prior work determined that efficient alphavirus replication requires catalytically active macrodomains; however, which form of the modification requires removal and from which protein(s) had not been determined. Here, we present evidence for the specific requirement of poly-ARH activity to ensure efficient productive infection and virus replication.

Comparative genomic analysis reveals varying levels of mammalian adaptation to coronavirus infections.

Severe acute respiratory coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is of zoonotic origin. Evolutionary analyses assessing whether coronaviruses similar to SARS-CoV-2 infected ancestral species of modern-day animal hosts could be useful in identifying additional reservoirs of potentially dangerous coronaviruses. We reasoned that if a clade of species has been repeatedly exposed to a virus, then their proteins relevant for viral entry may exhibit adaptations that affect host susceptibility or response. We perform comparative analyses across the mammalian phylogeny of angiotensin-converting enzyme 2 (ACE2), the cellular receptor for SARS-CoV-2, in order to uncover evidence for selection acting at its binding interface with the SARS-CoV-2 spike protein. We uncover that in rodents there is evidence for adaptive amino acid substitutions at positions comprising the ACE2-spike interaction interface, whereas the variation within ACE2 proteins in primates and some other mammalian clades is not consistent with evolutionary adaptations. We also analyze aminopeptidase N (APN), the receptor for the human coronavirus 229E, a virus that causes the common cold, and find evidence for adaptation in primates. Altogether, our results suggest that the rodent and primate lineages may have had ancient exposures to viruses similar to SARS-CoV-2 and HCoV-229E, respectively.

Identifying the Zoonotic Origin of SARS-CoV-2 by Modeling the Binding Affinity between the Spike Receptor-Binding Domain and Host ACE2.

Despite considerable research progress on SARS-CoV-2, the direct zoonotic origin (intermediate host) of the virus remains ambiguous. The most definitive approach to identify the intermediate host would be the detection of SARS-CoV-2-like coronaviruses in wild animals. However, due to the high number of animal species, it is not feasible to screen all the species in the laboratory. Given that binding to ACE2 proteins is the first step for the coronaviruses to invade host cells, we propose a computational pipeline to identify potential intermediate hosts of SARS-CoV-2 by modeling the binding affinity between the Spike receptor-binding domain (RBD) and host ACE2. Using this pipeline, we systematically examined 285 ACE2 variants from mammals, birds, fish, reptiles, and amphibians, and found that the binding energies calculated for the modeled Spike-RBD/ACE2 complex structures correlated closely with the effectiveness of animal infection as determined by multiple experimental data sets. Built on the optimized binding affinity cutoff, we suggest a set of 96 mammals, including 48 experimentally investigated ones, which are permissive to SARS-CoV-2, with candidates from primates, rodents, and carnivores at the highest risk of infection. Overall, this work not only suggests a limited range of potential intermediate SARS-CoV-2 hosts for further experimental investigation, but also, more importantly, it proposes a new structure-based approach to general zoonotic origin and susceptibility analyses that are critical for human infectious disease control and wildlife protection.